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GENERALIZATION OF BOUNDARY-LAYER MOMENTUM-INTEGRAL EQUATIONS TO 
THREE-DIMENSIONAL FLOWS INCLUDING THOSE OF ROTATING SYSTEM 1 

By Artur Maqer 


SUMMARY 

The Navier-Stokes equations of motion and the equation of 
continuity are transformed so as to apply to an orthogonal 
curvilinear coordinate system rotating with a uniform angular 
velocity about an arbitrary axis in space. ^4 usual simplifica- 
tion of these equations as consistent with the accepted boundary- 
layer theory and an integration of these equations through the 
boundary layer result in boundary -layer momentum-integral 
equations for three-dimensional flows that are applicable to 
either rotating or nonrotating fluid boundaries . 

These equations are simplified and an approximate solution 
in closed integral form is obtained for a generalized boundary- 
layer momentum-loss thickness and flow deflection at the wall 
in the turbulent case . 

^4 numerical evaluation of this solution carried out for data 
obtained in a curving nonrotating duct shows a fair quantita- 
tive agreement with the measured values. 

The form in which the equations are presented is readily 
adaptable to cases of steady , three-dimensional , incompressible 
boundary -layer flow like that over curved ducts or yawed wings; 
and it also may be used to describe the boundary -layer flow over 
various rotating surfaces , thus applying . to turbomachinery , 
propellers , and helicopter blades. 

INTRODUCTION 

The development of the boundary layer on the various 
parts of turbomachinery (compressors and turbines), heli- 
copter blades, propellers, and in curved ducts is influenced by 
centrifugal and Coriolis forces in addition to the pressure and 
viscous forces. As a result of these forces, the flow in the 
boundary layer not only has the characteristic velocity 
deficiency but also has, because of this velocity deficiency, 
direction different from that of the flow outside the boundary 
layer. Thus the behavior of the boundary layer in three- 
dimensional flow may be quite unlike the behavior in two- 
dimensional flow. The effect of these additional forces on the 
boundary layer has been realized for some time and the 
observed discrepancies in the boundary-layer behavior have 
usually been explained only in a qualitative manner as, for 
example, in references 1 to 4. 

The literature concerning the theoretical aspect of the 
three-dimensional boundary-layer flow is meager. For the 
laminar case most of the published work has been carried out 
in connection with the yawed wing (references 5 to 7). For 
the turbulent case, although a number of researchers have 
established the general form of the differential equations 


applicable, no actual solutions of these equations have been 
obtained. Tetervin, for example, presents boundary -layer 
momentum-integral equations in three dimensions for a fluid 
of variable density and viscosity (reference 8). Gruschwitz 
establishes the momentum -integral equations for boundary- 
layer flow along an arbitrarily curved streamline in reference 9. 
Burgers gives the differential equations on the develop- 
ment of boundary layers in the case of axially symmetric 
flows having a rotational component (reference 10) . Prandtl, 
in addition to presenting a form of three-dimensional 
momentum-integral equations, suggests the general procedure 
that could be followed to obtain a solution (reference 11). 
Experimental data are similarly lacking. In spite of con- 
siderable literature search, only the data of Gruschwitz 
(reference 9) for a curved duct and the data of Kuethe, 
McKee, and Curry (reference 12) for a yawed wing were 
found. 

As a result of research on this problem at the NACA 
Lewis laboratory, the boundary-layer momentum-integral 
equations are derived and presented herein for a set of 
orthogonal curvilinear coordinates, which may or may not 
be rotating about an arbitrary axis in space and can be laid 
out along a streamline of the potential flow. The so gen- 
eralized equations are then transformed by use of an assumed 
velocity distribution and friction law for turbulent boundary 
layer so that an approximate solution can be obtained for 
the boundary-layer momentum thickness and the direction 
of boundary-layer flow. Finally, a numerical solution is 
carried out for the Gruschwitz data in order to make a 
comparison between the estimated and actual measured 
values. 

The equations as given in their generalized form are 
readily adaptable to cases of steady, three-dimensional, 
incompressible boundary-layer flow involving centrifugal 
and Coriolis forces. The approximate solution, however, 
has been carried out only for the turbulent boundary layer, 
because in most of the aerodynamic configurations, where 
these equations apply, transition from laminar to turbulent 
flow occurs comparatively early in the flow process. A 
laminar form of the approximate solution can be obtained 
by simple substitution of a suitable velocity profile and 
friction law. 

It should be noted that whereas the differential equations 
describe the flow phenomena with only the accepted simpli- 
fications, the approximate solution depends to some extent 
on the assumed boundary-layer velocity profiles and the 
relation for friction. Both of these assumptions were made 


1 Supersedes NACA TN 2310, “Generalization of Boundary-Layer Momentum-Integral Equations to Three-Dimensional Flows Including Those of Rotating System” by Artur 
Mager, 1951. 

997810—52 1 



2 


REPORT 1067 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


on the basis of the data of Gruschwitz (reference 9) only, 
because the data of reference 12 were not adaptable to 
extensive computations for the purpose of this analysis. 
The measurements of Gruschwitz, on the other hand, have 
certain shortcomings as they were obtained in a nonrotating 
channel formed by two circular-arc shaped walls. Tlius the 
generality of the velocity profiles measured by Gruschwitz 
is in question. A revision of the approximate solution can 
therefore be expected when more data become available. In 
addition, any speculation on the occurrence of boundary- 
layer separation (which by definition is a special form of a 
velocity profile) would be absolutely meaningless; no further 
mention will therefore be made of this phase of the problem. 

SYMBOLS 

The following symbols are used in this report (the 
dimensions are given in right-hand column): 


A 

constant occurring in second ap- 
proximation for © 

(l- 1 ) 

a 

resultant acceleration vector in 
fixed (inertial) system 

(ft-*) 

a„ 

resultant acceleration vector in 
Cartesian coordinate system 

(it'*) 

B 

constant occurring in second ap- 
proximation for 9 

(0) 

I 

position vector of particle 

(D 

c 

curvature of x-axis (fig. 1), ^ 

r 1 ) 

d 

constant > 5 

/ 0.01255 P 1 dx\ 

(D 

E i(x) 

( 0} 

E-Ax) 

('*j P« "* A 

(0) 

E 3 (x) 

A 4 J*i u ) 

(0) 

?xx, ?xyt ■ 

rate-of-struin components 

(r l ) 

F 

resultant-force vector acting on 
particle 

(mlr 2 ) 

fxjrjz 

components of body forces per 
unit mass 

(IT*) 

G 

function describing boundary- 
layer velocity profile, also taken 

(0) 

g 

function describing boundary- 
layer velocity profile, also taken 

A' -IT 

(0) 


H,J,K,L 

quantities describing relations 
among various characteristic 

(0) 


loss thicknesses in boundary 
layer 


h\ } /&2j ^3 

transformation coefficients 

(0) 

l 

length 


MN 

parameter determining nature of 
boundary-layer equations 

(0) 

m 

mass 


P 

static pressure 

(ml-Hr*, 

PxxjVyy , • * 

. components of stress per unit 
area in Cartesian coordinate 
system 

(. ml-'r 2 ) 

1 0 

resultant velocity vector 

(ft' 1 ) 

R 

perpendicular distance of particle 
from axis of rotation 

(D 

Re 

d U 

Reynolds number based on 9 X) — — 

(0) 

r 

radius of circle 

(1) 

S 

total path length dx 

%) o 

(l) 

s 

arc length 

(l) 

t 

time 


u,v,w 

values of u , v, and w outside 
boundary layer 

(ft- 1 ) 

P OyV O 

velocities in Cartesian coordinate 
system 

(It-') 

u,v,w 

time averaged velocities in curvi- 
linear coordinate system 

(ft- 1 ) 

x ,Y,Z 

Cartesian coordinate system 

(l) 

*,y,z 

orthogonal curvilinear coordinate 
system 

(l) 


function used in transformation 

(l) 

a 

bound ary -layer deflection angle 
measured from direction of re- 
sultant skin-friction stress to 
direction of flow outside bound- 
ary layer 

(0) 

0 

angle between X-axis and tangent 
to x-axis 

(0) 

5 

boundary-layer thickness 

(l) 

5b K 

displacement thicknesses in three- 
dimensional boundary layer 

(l) 

€ 

m easu r e o f b ou n d a ry -layer deflec- 
tion, tan a 

(0) 


slope of characteristic line 

(0) 

e 

generalized boundary-layer 
momentum-loss thickness, 
6 x Re 1/4 

(l) 

0 X f 6 Z) &xit &ZX 

momentum-loss thicknesses in 
three-< 1 im ension al bou nd a ry 

layer 

(l) 
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X 

variable of function \p, 

(D 

V 

kinematic viscosity 

(«->) 

ii v> i 

components of vorticity vector 

(<-’) 

p 

density 

(ml 3 ) 

Gz, Try, . 

apparent stresses existing in tur- 
bulent flow 

(ml~'t~ 2 ) 

To 

shear stress at wall 

(ml~'t~ 2 ) 

4 

arbitrary function satisfying 
equation (39) and boundary 
conditions 

(j?^] r™™]) 


U) 

angular velocity 

(t~ l ) 

3 

£ 

3 

* 

3 

components of vector w in Car- 
tesian coordinate system 

(r 1 ) 

Wxt 

Subscripts: 

components of vector w in curvi- 
linear coordinate system 

(r 1 ) 

i 

initial value 


X 

2 -direction 


z 

2 -direction 


I, II 

order of approximations 



For Gruschwifcz data-point designations and streamline 



Figure 1 — Transformation from Cartesian coordinates X, Y, 7 to orthogonal curvilinear 
coordinates j, y, z. 


DERIVATION OF BOUNDARY-LAYER MOMENTUM-INTEGRAL 
EQUATIONS 

The equations for steady flow of a fluid having constant 
density are derived in a Cartesian coordinate system X y Y \Z 




(P) 

(a) Channel and measuring plate soon from below. 

(b) Measuring plate seen from below, showing point and streamline designations. Shaded 
sections Indicate regions of potential-flow breakdown. 

Figure 2. — ExixTimental setup of nnischwitz (from figs, i ami 5 of reference !1). 

rotating with uniform angular velocity about an arbitrary 
axis in space. These equations are then transformed to an 
orthogonal curvilinear coordinate system x,y,z such that 
the x-axis can be placed along any convenient path in the 
XZ-plane, which is considered as a plane of a wall. These 
equations are then simplified in a manner consistent with 
the boundary-layer theory. If the path is chosen so as to 
match a streamline of the potential flow, only one velocity 
will exist outside the boundary layer, that along the stream- 
line. Furthermore, the changes in boundary-layer quantities 
in a direction other than that along the streamline are ex- 
pected to be relatively small in comparison with the changes 
along the streamline. Additional simplifications may thus 
be possible. Finally, integration through the boundary 
layer gives the generalized form of momentum-integral 
equations for three-dimensional flows that may or may not 
involve rotation of the system. 

Equations for steady flow of fluid with constant density 
in rotating Cartesian coordinate system. The Wavier- 
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> r Vv i 




• r ;^T.7;i 


X±A 

;^;A 


Stokes equations of flow for a fixed Cartesian coordinate 
system X,Y,7 (reference 13, p. 576) are 


DU. 
Dt 

DV. 


r , dpxx I &Prx i upzx 

=pfx+ -bX' + -dY' + - 


bpz 
' dZ 


Dt 

Dir, 

Dt 


4 , bp X y , dp ry , dpzr 

dZ + dF + dZ 

— V i dpxz i bprz . bp zz 

~ pJz+ 'bX + dF + bZ 


and the equation of continuity is 

bu 0 , bv a , dm 


dZ 1 dF 1 dZ 


= 0 


(la) 

(lb) 

(lc) 

(2) 


. •,,* ’..1 
nan'll 

- 

n *"• » 
X*X\-V'3 


V.V ‘- *-< 


It is now assumed that this Cartesian coordinate system is 
rotating with a uniform angular velocity a> and that the 
observations of the motion of the fluid particles are still 
made from a position rigidly attached to the same system. 
The velocity q 0 and acceleration a 0 are as seen by the ob- 
server, that is, they are relative to the X,Y,7 system. 
Because of the rotation, however, the X,Y,7 system is not 
an inertial system (reference 14, p. 53) and thus the second 
law of motion holds only with respect to acceleration a 
relative to some other system that is nonrotating, 

ma = F 

In terms of a 0 then (reference 14, p. 104), 

ma 0 + m,uX (o>Xb) + 2mwXq 0 — F 

Here o?a>X(ajX6) represents the centrifugal force and 
2m<j)'Xq 0 is the Coriolis force. 

Thus for a Cartesian coordinate system rotating with a 

DU 0 DV 

uniform angular velocity w, the expressions for — ^ 
DW 0 

and Y) Y mus ^ be modified by proper components of the 

Coriolis and centrifugal accelerations. For steady flow, the 
component accelerations as referred to a rotating Cartesian 
coordinate system are therefore 

DTJo JJ 2) H o | T7 ^ C o , TJ7 t) Uo [ o/.. TJ 7 T7\ 

— C o T V 0 ^ yr B o + 2(a>y rr 0 <j)z Vo) W /? 

(3a) 

Dt ~ Lo dZ + K 57 J rW 0 ^g+-{<o a L'— wx-W,)— « R^y 

(3b) 

DW. d W. dTF„ w dTF, bR 

Dt ~ L ° dZ + dF + ’ dZ w r L.)—u >‘R^ 

(3c) 


dZ 


dF 


The equation of continuity, which does not involve any 
accelerations, remains the same. 


Transformation to orthogonal curvilinear coordinate sys- 
tem, — Transformations similar to those of Gruschwitz 
(reference 9) are used as indicated in figure 1 with the pre- 
caution that the system remain right-handed. - 


r =J cos 0 dx + z sin 0 
Z=Z l + 2 cos 0 


(4) 


where 


and 


Z v = constant- 


■JT 


sin 0 dx 


0 = 0(x) 


Use of these transformations permits an arbitrary curva- 
ture of the x-axis in only one plane, the -XT-plane. Thus 
the solution is somewhat restricted. In two-dimensional 
boundary-layer investigations, however, it is found that the 
boundary-layer equations are unaffected if the radius of 
curvature in the -XT-plane is large as compared with the 
boundary-layer thickness (reference 15, p. 120). In three- 
dimensional boundary layer the same limitation will 
probably apply providing, of course, the values of oj x , 
and o) z are properly adjusted to take care of this additional 
curvature. Setting 

c= ~^ (curvature of x-axis) 


gives 


^=(l+C2) COS 0 


dF 

dx 

dZ 


= 0 


dZ 

by 

bY 

by 


= 1 


dZ 

bz 

bY 

Oz 


— sin 0 
= 0 


^ =— (1 -\-cz) sin 0 ^-^=0 ^?=cos 0 
ox Oy Oz 

The elements of length at (x, y, z) in the direction of the 
increasing coordinates are (reference 15, p. 101): 

hidx, h 2 dy , hzdz 

Thus, 

(dsy= (k,) 2 (dx) i +(h 2 y(d y y+ (h 3 y 2 (dz) 2 = (dX)*Y(dYy+(dZ)- 2 

But because 

dX — — dy + ^- dz 

Ox Oy 9 Oz 


and so forth, 


and 


(ds) 2 =(l+czY (dx) 2 +(dy) 2 +(dz) 2 

— ( 1 — {— C3) h>= 1 A 3 = 1 


05 r 



GENERALIZATION OF BOUND ARY-LAYER MOMENTUM-INTEGRAL EQUATIONS 

1 


The expressions for the linear accelerations can be written 
directly, as given in reference 13 (p. 158 ). (It should be 
noted that the h values herein are reciprocals of those in 
reference 13 .) The components of a gradient now are 

_L JV ' l_b_ 
h 2 by h z bz 


1 A 

hi bx 


whereas the components of remain 

V (j) t — W toy 11 Wj — L <J&z Tj U>y — Vto x 


Thus the accelerations in the rotating x t y,z 
written as 


system are 


DU_ U bU V bU W bU 

h\ bx h 2 by A3 bz 


Dt 


~r(u tt 1- v 7 ^)+ 

,A 2 \ dy bx J 


W 

A 1 A3 




bR 

bx 


And the expressions for and follow from symmetry. 


The equation for the divergence now has the form 
1 


div q 0 z 


A \ h 2 A3 


[A (A 2 A 3 L0+^ (MitO+sj (A.A.W) 




whereas the components of the curl q 0 are 

i - JsK W > 5 ' W ] 




In order to obtain the viscous terms the preceding expres- 
sions are used in the expansion of 


v [grad (div q 0 ) — curl (curl q 0 )\ 


If equations ( 5 ) are substituted into these general expres- 
sions and the differentiations are carried out, the equations 
for flow in an orthogonal curvilinear coordinate system rotat- 
ing with an angular velocity w are obtained. The body forces 
are neglected here. 


1 +zc bx' by bz 


1 +zc 


UW - 


1 


1 +zc bx 


2 (w v IF — oitV)"- 


1 1 bP 


bU dc 


1 -\~cz p 

b 2 U , b 2 U 


bP , r 
5 ? + l( 


(I+C2) 3 bx dx" 1 " dy 2 " 1- bz 2 "*"l +cz bz (1 +cx) 


bU Uc 2 


W dc , 2c bWl 


(1 +C 3 ) 3 dx ' (1 +czY bx J 


(6a) 


utr+rg+wg- 

bx by bz 


——JR %£+2 («,U-» z W) 

1 +zc by 


b 2 V 


l&p , r i 

— p dy i ’”L( 1 + cz ) :! dx2 U+c 2 ) 3 a x dx 


bV dc 


b 2 V . b 2 V 


by 2 

U 


bz 2 1 

dW 


1 +cz bz 


I ] 


(6b) 


■V^.+ W iW 


bz 


1 +zc bx n ' by 

2 <„,F-*0) ;|£+«[; 

U dc z bW dc 


Id - zc 
1 


C*— ’*§ ?+ 


b 2 JV 

(1+C2) 2 bx 1 


(1 +cz ) 3 dx (1+C2) 3 dx dx 


&w d-ir 

dy 2 d2 2 


Wc 2 


(1+C2) 


-«+ 


bW 


m 


1 +cz bz 
1 


(1 +cz ) 2 bxj 


bU + bV 


bW , Wc 


bz 1 +cz 


= 0 


( 6 c) 

(6d) 


1 +cz bx by 

In the general orthogonal coordinates, the expressions for 
the rate-of-strain components are 


(± bU . V bht JF_ bkC\ 

bx^hiht by^IJh bz J 


by 

W bh 2 


• =2 ( 


■2 dy 
i diF 


. U dAA 

A 2 A 2 d2 AjA 2 dx / 


U bhj V bh 3 \ 


h 3 bz t A 5 A, dx A 2 A 3 dy/ 

= A3 A A 

A 2 dy \ A3 / A3 d 2 \A 2 / 


A, A 
A 3 d 2 


/£\ , A, A^ 

\Ai/ '"hi bx\h 3 J 

„ _^A/YAj-ZJL/TA 

A! dx \A 2 /^" A 2 dy \A|/ 


The viscous terms in equations (6a), (6b), and (6c) may be 
expressed using the rate-of-strain components as 

r l b 2 U z bU dc , b 2 U , b 2 U 

v Lc 


J1+C2) 2 dx 2 (1+C2) 3 dx dx^by 2 ^ dx 2 + 


bU Uc 2 


W dc . 2c d UH 


1 +C 2 dx (l + cx^Al+cxfdxAl+cx ) 2 dx J 


1 b 2 U 

r 1 de„ 

^ € vx j 

, d«« . 

2 ce zl " 

^(1 + cz ) 2 bx 2 

^Ll+C2 bx 

h dZ H 

h d 7 ^ 

1 d" C2_ 


■[ 


d 2 F 


bVdc , d 2 F , d 2 F 


,( 1 +C 2 ) 2 dx 2 (l+cx) 3 dx dx 


{ 


be 


VV [ 

by ^ J +c 


1 be, v , be 


bx ^ bz ^1 +cz 


h- 


by* 


dz‘ 


( 7 a) 

bVl 


1 +cz bz J' 
( 7 b) 



6 


REPORT 1007 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


’[(1+c 


b 2 W 


U dc 


bWdc 


_(1 -\-czf bx 2 (1 +cz) 3 dx (l+cz) 3 dx dx 
l ) 2 W , b 2 W Wc 2 , c bW 2c 


by 1 

= "[tT 


bz 2 (1+C 2 ) 2+ 1+C2 bz (1 + cz) 2 


be r 


bc yz i b€ z , 


cz bx by bz 1 +cz 


(&ZZ &zz) 


I 


(7 c) 


Equations (6) are directly applicable to the laminar flow. 
For turbulent flow, because of the velocity fluctuations it is 
necessary to modify the stresses by addition of the so-called 
Reynolds 7 stresses. Thus, making use of the parallel form 
in equations (7), the Navier-Stokes equations of motion for 
turbulent flow may be written in terms of the apparent 
stresses as 


^i+ v *E + w 

l+cz bx T by ^ bz^l+cz 
1 u 2 R^-+2(u v W-u,V)-- 


rotating orthogonal curvilinear coordinate system evaluated 
at 2 = 0 , 

bu bu bu IbP 2T? bR , / b 2 u\ 

u ^ +v ^ +w ^ +cuw= -^ +uR ^- 2 ^ w+ \W 2 ) 

(9a) 
(9b) 


2E> bR , of \ 1 bP 

— <ji a zr — 2 (oi z u — w) = — 

by p by 


bw , bw . bw , 1 bP , ^bR , n 

u te +v ^ +w ^- uc= --pte + "- R ^+ 2 “' u ^ v KwJ 


fbhv\ 
by 2 ) 
(9c) 


UW — 


1 


1+c 

IU 

pLi + 


dx 

d<r ; 


I . 
. ^ „ 1 


1 + C2 p dx 




dx by bz 


2cr„ ~ 
1 + C2_ 


(8a) 


TT bV _ .bV , TT-dF 1 jptlff , T r TT r. 

U -r — (- V — (- H -r — r — u‘R r — r 2 (ojj u — bi z W ) 

dx by bz 1 + cz by 


= 1 dP 1 

p by p 


w 


drjy . b<J v 


cz dx by bz 


| CT yt "1 

d2 ' I + C 2 J 


(3b) 


for the laminar case. For the turbulent boundary layer, a 
corresponding set of equations is obtained with the substitu- 
tion of - for v ( and — ^-^for v 

p by \by 2 / p by \by 2 J 

Equation (9b) shows, as pointed out in reference 10, that 
because all the terms on the left-hand side of the equations 
are of the order of magnitude of one, within the boundary 
layer, P can vary at most by an amount of the order of 5. 
It is reasonable then to neglect this variation and consider 
P solely a function of the flow outside the boundary layer. 
Thus, if x is chosen to coincide with a streamline of the flow 
outside the boundary layer, V=W= 0, and by integration 
of equation (9a) with the effect of viscosity neglected the 
following relation is obtained: 


P=constant— - pU 2 +^ P^R 2 


( 10 ) 


U bW . „bW . ytrbW 


c 

l+cz 


U 2 -u 2 R^+ 

d2 


2 {wxV — U ) = 


1 bP ,1 [" 1 br zz 

p dz ^ p Ll +cz dx 


which is a form of the equation of Bernoulli. 

Furthermore, because outside the boundary layer the flow 
with respect to some nonrotating set of coordinates is iiTOta- 
tional with reference to the rotating coordinates the com- 
ponents of the vorticity vector become 


br y2 . b<r z 
by bz 



(3c) 


Simplification for flow within boundary layer. — Equa- 
tions (6) and (8) are equivalent to the complete Navier- 
Stokes equations. Within the boundary layer, however, 
certain terms whose contribution is relatively unimportant 
can be neglected. If the y-axis is taken as normal to the 
wall, the boundary-layer flow then takes place over the X2- 
plane (or the ATZ-plane). All terms are now made dimen- 
sionless by referring the lengths to some body length, the 
velocities to their free-stream values, and so forth, as ex- 
plained in reference 16 (p. 45), and all quantities of the order 
of magnitude of <5 or smaller are neglected. Furthermore, 
because the boundary-layer flow along a definite path 
2=0 is of interest, additional simplifications are possible. 
Setting 2=0 restricts the equations, because the general 
boundary conditions (not on the x-axis) cannot be satisfied. 
It will subsequently be seen, however, that these general 
boundary conditions are unnecessary in the solution of the 
final equations. These simplifications yield the Navier- 
Stokes equations for flow within the boundary layer in a 


{ = — 2oj 7 "1 

V = L (II) 

r = — 2 J 


This assumption of irrotationality is not always true and in 
some applications, such as the later stages of an axial com- 
pressor, it cannot be used. As long as vorticity is distributed 
according to some definite pattern, however, a relation be- 
tween the components of vorticity and the components of 
rotational velocity may be found and substituted for equa- 
tions (11). 

Substituting again in the expression for the components 
of vorticity gives 


And for 2 = 0 , IF=0, which is along the streamline, the 
expression for curvature becomes 


c 


/2^ , J _b£\ 
bz) 


( 12 ) 



GENERALIZATION OF BOUNDARY-LAYER MOMENTUM-INTEGRAL EQUATIONS 
The equation of continuity remains 

(13) 


^ 0 

dx oy oz 


Generalized boundary-layer momentum-integral equa- 
tions.— In order to obtain the boundary-layer momentum- 
integral equations, equations (9a) and (9c) are integrated 
with respect to y through the boundary layer to some con- 
stant height d such that 


d>& 


dz 


j* u'dy + '^f'wdy+^f'w dy-V §- r j*udy— 

u tzl dwdy+ hS! uwdy ~ 4 v$! uwdy - 

2 b ikS* uwdy+2o> “.i wd y =Ud ^~ T ~f (14a) 


and 


g-j‘„viy+ *-f‘„-dy+c J t Wd y-c^ uhly 


(14b) 

Those equations apply equally well for the laminar oi tur- 
bulent boundary layer, with the value of r 0 representing the 
shear stress at the wall accordingly adjusted. By suitable 
use of equations (12) and (13), these equations may be 
transformed to 

V Iff' JV— Xxiv-A-’ »»*> 


and 


(15b) 

The following definitions are now introduced: The mo- 
mentum thickness in the z-direction of the flow in the x- 
direction, 

e z= ^j\u-u)udy (16) 

The displacement thickness in the z-direction, 

K=jjj\u-u)dy 


(17) 


The momentum thickness in the z-direction of the flow in 
the z-direction, 


'■“ 77 * t vMy 

The displacement thickness in the z-direction, 

1 C d 

d ’=Ujo Wdy 


7 

(18) 

(19) 


The momentum thickness in the z-direction oi the flow in 
the x-direction, 

( 20 ) 


1 C d 

9 xz =jj- 2 I ( U—u)wdy 


The momentum thickness in the x-direction of the flow in 
the z-direction, 

% d 

(21) 


wudy 


u,x ~U i 

All these thicknesses, as in two-dimensional boundary- 
layer theory, have a dimension of length. Furthermore, 

{U-u)wdy=^ l ^ uwdy=e tI 

( 22 ) 

With the use of definitions (16) to (21) and equation (22), 
equations (15a) and (15c) reduce, for z=0, to 


and 


<)0 Z . 5(6 a d xz ) . 1 &U , n g*\ | 2 dL f ^ ^ \ — 

^1+ dx“ + t/dT {dl ~ 91 5 ' )+ U bx {5 ‘ 9 " } 


^(0r+0d = 


' pU 2 


(23 b) 


Reduction of equations to forms obtained by other 
investigators. — If only two-dimensional flow exists, that 
is, if c = 0, «?=0, and w=0, then equation (23b) vanishes and 
equation (23a) becomes an ordinary Karm&n momentum- 
integral equation 

dftr I 1 ^LL(Of) 1 sq T °» z - 

If w=0, that is, if the system is nonrotating, equations 
(23a) and (23b) become identical with the equations of 
Gruschwitz (reference 9). 

Setting c = 0 in equations (9a) and (9c) makes these equa- 
tions identical with the equations of Burgers (reference 10), 
who carried out his derivation for a Cartesian coordinate 
system. 

Finally, if the system of coordinates is chosen so as to 
maintain the right-hand rule and c is set equal to —* thus 

establishing the x-axis as a circle, then^ = l and because 

of axial symmetry all derivatives with respect to x vanish. 
The coordinates are now assumed to be in a fluid that is 
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motionless at great distance from the surface of the rotating 
immersed disk. Thus, 

U=W=0 w =0 

Integration of equations (9a) and (9c) gives, after some 
manipulation, 

2 *i( r '-$l uwd y)=- T -j l2 « ri 

and 

which are identical with equations of Karm&n for the 
rotating disk (reference 17). 

APPROXIMATE SOLUTION OF MOMENTUM-INTEGRAL 
EQUATIONS FOR TURBULENT BOUNDARY LAYER 

Transformation and reduction of dependent variables. — 

In order to obtain a solution of the momentum-integral 
equations, additional relations are needed describing the 
velocity profiles existing in the boundary layer and the 
friction at the wall. 

With the use of a suggestion by Prandtl (reference 11), 
the expressions for u and w that will be used are 



with boundary conditions on 0 and g 


for moderately large e. In addition, in order to check this 

assumption, the value of — ™ for several experimental 

« U 

velocity profiles and values of e ranging from 0.216 to 
0.670, as obtained from reference 9, is plotted against y/5 in 
figure 3. The results of this plot indicate indeed that Gg is 
independent of «. 
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Figure 3.— Plot of ! ^ against ! for various data from reference 9. 

< U i 


In parallel to the two-dimensional boundary-layer theory, 
the following definitions are made: 


for y= 5 , 

C=l, 

5=0 


for y= 0 , 

cs 

If 

© 

5 = 1 


and with < defined as 

€ = tan a 


(25) 


where a is the angle between the direction of the resultant 
skin-friction stress and the direction of the flow outside the 
boundary layer. Because of this definition of e, g=l at 
y = 0 because 

dw 

dy. 


i/^o oil 

dy 


or 


(26) 


Mathematically, such use of t implies a linear variation of 
w with e and makes possible the dissociation of the w velocity 
profile from its scale and direction. Because the flow must 
change direction in the boundary layer from that at the wall 
to that in the free stream, there is no reason to assume that 
such a dissociation is actually possible. In other words, 
there is no reason to believe that g should be a function of 
(y/8) alone and not of e as well. In accordance with reference 11, 
however, this approximation is certainly admissible for 
small values of < and gives results of qualitative accuracy 


f (1— G)dy 

JO _ TT 

Pd 

J (I —G)Gdy 


(l-G)Ggdy 


L 

(l—G)Gdy 

f'Ggdy 
f\l-G)Gdy 

f*GYdy 
J\l -G)Gdy 




K 


(27) 


The relations among the various thicknesses may then be 
written 


Sjz == € jBx 

5 ] = tKe z 

e,=t 2 L9 XJ 


( 28 ) 
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GENERALIZATION OF BOUNDARY-LAYER MOMENTUM-INTEGRAL EQUATIONS 


The quantities ff, J, if, and L are functions of G and g. 
Because G and g are representative of the velocity profiles in 
the boundary layer, the changes in these velocity profiles 
must be reflected in turn in the values of H, J, if, and L . 
In other words, the external forces acting on the boundary 
layer and influencing the changes in the shape of the velocity 
profiles also cause a variation in H, J, if, and L. Unfortu- 
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(a) Experimental velocity profile from reference 9 (point 15, III). Data obtained in curved 

duct; 40 millimeters. 

(b) Experimental velocity profile of reference 12 converted tox, y, z coordinate system. Data 

obtained in boundary layer of yawed wing. 

Fio.ube 4.— Comparison of assumed O and Og with experimental velocity profiles. Assumed 
profiles: <3— (v/«*) VT ; Gi~{yll)W (1 -?/$)*. 
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Fic.ube 5. — Comparison of assumed correction function g(yl&) with data of reference 12 con- 
verted tox, y, z coordinate system. Assumption: f(y/J)*( 


nately, the avaflable data of reference 9 do not involve 
large" changes in the shape of the velocity profiles and the 
quantities H } J, if, and L. This behavior of the velocity 
profiles is verified in figure 3. The data of reference 12 do 
indicate large changes in the shape of the velocity profile, 
however, the data are not presented with sufficient detail to 
permit an accurate evaluation of H , J, if, and L. Thus, 
until more extensive experimental data become available, 
the quantities H, J } if, and L are. assumed to be constants 
that can be evaluated either by assuming a suitable form for 
G and g or by computing directly from Gruschwitz data. 

In accordance with reference 9, good assumptions for G 
and g are: 


a-l 




(29) 


An indication of the degree of fit afforded by these expres- 
sions can be obtained from figure 4 (a), where a calculated 
profile with n=7 is compared with one of the profiles of 
Gruschwitz. Other profiles of Gruschwitz data give similar 
results. It should be noted that this good agreement should 
not be interpreted as meaning that assumptions (29) will 
always give a good representation of the velocity profiles in 
the three-dimensional turbulent boundary layer. Figure 4(b) 
shows a comparison similar to that of figure 4(a) with pro- 
files converted to the x,y,z system using data from reference 
12. Equations (29) do not afford a good fit in figure 4(b), 
although the equations do represent the general behavior 
of the velocities. This comparison is further illustrated in 
figure 5, where the value of g (yl&) as obtained by converting 
the profiles of reference 12 to the x,y,z system at indicated 
points is compared with (1 —y/S) 2 . 

With the use of relations (29), H, J, K, and L are com- 
puted as 


H= 


2 + n 


71 


J= 

K- 


n, 2 (l ln-j-7) 


L= 


‘(2n+l)(3n+l)(3n+2) 
27l 2 (2 + 7l) 

\'2n+l)(3n+l) 

6 a 1 


(3n+2)(2»+l)(5n+2)J 


(30) 


which for n=7 give 


H= 1.2857 
J= 0.5423 
K= 2.6727 
L= 1.1285 


(30a) 
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Averaging the values along line III of Gruschwitz data 
(fig. 2) results in 


H= 1.37 " 
J= 0.550 
K=2A3 
L= 0.968 , 


(30b) 


This relatively good agreement between the two sets of 
values is also indicative of the over-all fit of the assumed 
expressions for G and g to the data of reference 9. 

The additional relation that is needed for the solution of 
the momentum-integral equation is the expression for surface 
friction. In reference 9, Gruschwitz demonstrates that 
Karmen's friction law 

Up- 001255 (A)* (31) 

appears to be valid in the three-dimensional boundary layer 
as well. Substituting relations (26), (28), and (31) in 
equations (23) yields 


T+ J9 ’M- 4 v Jte ’ 




As shown in the appendix, these equations can be either 
hyperbolic, parabolic, or elliptic, depending on the shape of 
the velocity profiles existing in the boundary layer. For 
u=U(y/S ) 11 * and g—(l— y/5) 2 , the equations are always 
elliptic. 

Simplification of equations and approximate solution. — 

The relative importance of the various derivatives in equa- 
tions (34) is now determined. First, 9 and * are assumed to 
be quantities that are smaller than one, which can be accom- 
plished simply by referring all lengths to a total path length 
S and a to 45°. As a result of this assumption, all deriva- 
tives in 0 and « become of the order of magnitude of 0 or t. 
Equation (34a) is then divided by 4/5, so that the coefficient 
of bB/bx is 1 and the coefficients of bBjbz and btfbz are Je 
5 

and ^ JB, that is, of the order of magnitude of e and 9, 

respectively. In a similar way equation (34b) is divided by 
(K— J)B so that the coefficients of d9/c)z and bej’bz become 

4 

^ e/0 and 1 (order of magnitude of one), respectively, and 

4 € 2 

those of bB/bz and btjbz become t L ~ and 2 L t (order of 

o y 

magnitude of c), respectively. Then, if e is small as com- 
pared with tan 45° and 9 is small when compared with S, all 
terms of the order of magnitude of € 2 , 9e, and 9 2 may be 
neglected, which gives 


and 

2L >* h+ L '' | 57+ 


pU 2 

(32b) 


be 

bx 


+ 


r/5f?+9\ 

1 bU~ 

Lv 4 ) 

1 U bx_ 


.01569 


(35a) 


and 


,9 i ar 0.01255 n 
+ 5 V bx + (K- J) 0j * 


b^rnbe 

dx |_5 0 bx 


(35b) 


Because of the form of the relation for friction, an advan- 
tageous transformation of variable is 

B=d z Re l/i (33) 


in order to eliminate the Reynolds number from the equa- 
tions. 

With the use of equations (33) and (12), two nonlinear 
partial differential equations for 0 and c applying along z— 0 
are obtained from equations (32), 


4» + ^«“ +J e|i+ 

o bx b bz bz 


1 d?7 . 
bx^ 


iL y-r-Ju 


Y Ct — J* y?"] 9 — 0.0 1255) = 0 (34a) 

5 5 L J ) 

and 


I (/£• — D« H+| £<’ H+(K- 57+2M- §;+ 


(^« 2 -f)^ 0 + 0.01255fJ=0 (34b) 


These two expressions show that the primary changes in 
0 and e in equations (34a) and (34b) occur only in the 
z-direction and thus the description of the phenomena only 
at z= 0 is justified. 

A solution of equations (35) can now be obtained by suc- 
cessive approximations because U, uy, and c are assumed to 
be known functions of x. First, equation (35a) is solved, 


0i=0, 


/ 5 W+ 9 \ 

/£7y 4 } Q.Q1569 r x 

{u) 


/5H+9\ 

CA 4 ' dx 


(36) 


The values of 0i(a) are then used in an approximate solution 
of equation (35b) 


/e ( WDiY i . 

e,_<< W \u) E t (x) + 


(J-K)(QiyU* E^x) 


where 


. JT(1 + H) c +H (0,/ U* E x (x) dx 

0.0125 5 rx dx 


(37) 

(38) 



GENERALIZATION OF BOUNDARY-LAYER MOMENTUM-INTEGRAL EQUATIONS 


With €j(x) known, the approximation for 0 could in turn 
be improved by again solving equation (35a). If a grows 
large along the path, however, it is more advantageous to 
consider the following equation : 

|+^|f+[( ? ^)r^+I'-T' 7 -F] 9 -o 01569 

(39) 


It is thus hoped that neglectionof befdz will not affect the ac- 
curacy of the solution to a very marked extent. The solution 
of equation (39), which may be obtained by the method of 
Lagrange, is 


where 


and 


e,I= [*( 2- ' 7 k << * X ) 


E 3 (x) 


/5M+9N 

•v « > E,(x) 


+ 


/ 5//+9 \ 

0.01 569 E 3 {x) f r tA 4 > Ej{x) 
/SfffSN J x E 3 (x) 


dx 


Ctdx 


-f r 

E 2 (x) = e i ^ 


IS r f r , 

yJ < 77 ^ 
Et(x) = e* L 


(40) 

(41) 

(42) 


is an arbitrary function satisfying the boundary condition; 
when j — Xi and 3=0, then 0 n = Oi. Setting 


(^z — J j* t (1x^=1 


at x=-Xi gives X=0. In addition, for z — 0, X is of the order of 
magnitude of e. 

Expansion of $ in Maclaurin’s series about X(x<) yields 


^(X)=^(0) + X^'(0)-)-^j ^(0)+ • * • 


Because there is only one boundary condition, it is possible 
to determine only one of the constants in this expansion; 
consequently, $ cannot be uniquely established. The fact 
that X is of the order of magnitude of e, however, suggests 
that the assumptions made for ^"(0), *'"(0), and so forth, 
are successively less important. Thus these derivatives may 
arbitrarily be expressed by a single constant, 


/5W9\ 

*(X)=0,E7A 4 ^1+^4X + 




/ 5gj-9 \ 

U\ 4 ) e AX 

(43) 


where A from purely dimensional considerations must have 
the dimensions of l~ l . From expressions (41) and (42), it is 
suspected that 

(44) 




where B must be obtained from the experimental measure- 
ments. 

It should be noted that because ^(X) cannot be uniquely 
determined other functions of X satisfying the single boundary 
condition could be used as well. The function e xx is chosen 
only because it is convenient to use and parallels the expres- 
sions (41) and (42). This arbitrariness of the functional form 
of and the value of .4 is due to the consequences of assuming 
2 = 0 , and thus it is probably not advisable to carry any 
further approximation for «n and so forth. 

In solving equations (36) to (44), either set of values for 
H , K y and L may be used. Because the averaged values 

(30b) are probably more representative, having been obtained 
by evaluating experimental data at a number of different 
positions, it is advantageous to use these values in com- 
putations. 

COMPARISON WITH EXPERIMENT 


In order to check the validity of the approximate solution, 
the boundary layer along four streamlines of reference 9 was 
computed and compared with the measured values. The 
designation of the streamlines and data points is illustrated 
in figure 2. Because the data were taken along curves I to 
V of figure 2, the computation along a streamline requires 
first an interpolation among the various data points. As a 
result of this interpolation, the computations could not be 
carried through the full length of each streamline. X alues 
(30b) were used for quantities H , J, K, and L. The constant 
B was obtained by fitting along streamline B the solution 
for 0 n , so that at x=S, 0 n ~0 measured. In this manner, 
the value of B was found to be 38.5. This value was then 
used in computations of streamlines A, C, and D. It is 
noted that B=^S.5^7(Re l/i ) i} although justification for such 
a dependence cannot be made. In all integrations Simpson s 
rule was used. 


The results of the computations arc plotted in a nondi- 
mensional form and compared with the interpolated measured 
values in figures 6 and 7. A study of these figures reveals 
a fair quantitative agreement between the measured and 
estimated values of 9 and a. As the values of e— >tan 45° 
(fig. 7) the first approximation for 0 in figure 6 becomes 
progressively worse, which is remedied by the second approx- 
imation. The poorest agreement is obtained along stream- 
lines A and D, which because they are closest to the walls 
might be affected by the flow in the corners of the duct. 
Streamline D especially may be affected inasmuch as 
Gruschwitz mentions the existence of separation on the 
convex wall. 

The fair quantitative agreement with the measured values 
is not to be interpreted as a conclusive check of the validity 
of the procedure and the assumed values in all cases of three- 
dimensional boundary-layer flow. The suggested procedure 
simply represents the best that can be done in view of the 
meagemess of the available data. Because the Gruschwitz 
data do not involve the effects of uniform angular velocity 

and because the variations in — and jj are small, 
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(a) Streamline A; 5—22.8 inches; 9i«0.292 inch. 

(b) Streamline B; 5-34.75 inches; 9,-0.289 inch. 

(c) Streamline C; 5-37.8 inches; 9,-0.286 inch. 

(d) Streamline D; 5=21.25 inches; 9,-0.244 inch. 

Figure 6.—Comparison of calculated and measured generalized momentum thlcfcness. Experimental data from reference 9. 
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it could be maintained that this check of the procedure has 
been carried out on a somewhat special case. For that 
reason, it is desirable that additional experiments be carried 
out in setups that eliminate the present shortcomings. A 
larger variation of Reynolds number should also be used. 
With additional experiments, a modification of the values 
of H , J, K, L , and B, together perhaps with some refinements 
of the procedure, will be in order. It might be well to remem- 
ber, at such time, that because of the necessary empiricism 
involved (which results from the very limited knowledge of 
turbulent phenomena), long and tedious computations would 
rarely be worthwhile. 

CONCLUSIONS 

The following conclusions can be drawn from an analysis 
of the three-dimensional momentum-integral equations and 
a comparison of the numerical results with the Gruschwitz 
data for turbulent boundary layer: 

1. Within the boundary layer the static pressure can vary 



at most by an amount of the order of magnitude of the 
boundary-layer thickness <5. 

2. It is possible to generalize the velocities in the boundary 
layer by use of two characteristic quantities 5 and e where 
t is the tangent of the angle enclosed by the direction of the 
resultant skin-friction stress and the direction of the flow 
outside the boundary layer. 

3. When the generalized boundary-layer momentum-loss 
thickness 0 is small as compared with the total path length 
and € is small as compared with tan 45°, the primary changes 



(b) Streamline B; £-34.75 inches. 


(a) Streamline A; S— 22.8 inches. 

Fioube 7.— Comparison of calculated and measured boundary-layer deflection at wall. Experimental data from reference 9. 
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in 0 and t occur along the streamline of the flow outside the 
boundary layer. 

4. The three-dimensional boundary-layer momentum- 
integral equations can be either hyperbolic, parabolic, or 
elliptic, depending on the relative magnitude of the parameter 
MN, which in turn depends on the shape of the velocity 



profiles existing in the boundary layer. The power-law 
profile when used with the correction function g = (l—y/&) 2 
always results in elliptic equations. 

5. The approximate solution of the three-dimensional 
momentum-integral equations shows a fair quantitative 
agreement with the values measured by Gruschwitz. 

6. Additional experimental data are necessary to establish 
more generally applicable values for form parameters H, 
J, K, and L and B , the constant used in the second approxi- 
mation for 0. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, November T } 1950 



(c) Streamline C; 5*37.8 inches, (d) Streamline D ; 5*21.25 inches. 

Figure 7. Concluded.— Comparison of calculated and measured boundary-layer deflection at wall. Experimental data from reference 9. 


APPENDIX 


ADDITIONAL REMARKS ON THREE-DIMENSIONAL 

BOUNDARY-LAYER MOMENTUM-INTEGRAL EQUATIONS 

In order to obtain the approximate solution of equation (34), 
it was shown by comparing the relative order of mag- 
nitude of the coefficients that some of the terms may be 
neglected. Care must be taken with such simplifications 
inasmuch as various implications of the equations in question 
may be obscured by this procedure. For this reason, aside 
from the approximate solution, the character of equations (34) 
was also investigated in detail. 

By use of the procedure outlined in reference 18 (p. 38), 
along 2=0 the system of equations (34) is found to be 
hyperbolic when MN^> 1, elliptic when MN<^ 1, and para- 
bolic when MN= 1, where 


MN as 




t K-J)J e t ,e. 


(45) 


Because J, K , and L are functions of G and g, the character 
of equation (34) depends on the shape of the velocity profiles 
in the boundary layer. 

It should be noted that when MN=0> then L= 0, which is 
only possible if #=0 or £=0, and in turn u = 0 or w=0. If 
the trivial case u— 0 is neglected, it is established that when 
w=0, € = 0 as well. But for e = 0 and w= 0, equations (34) 
reduce to a special case 


and 


4 50 e bU 

5 dr U bx 


( 5g 5 ~t 9 )= 0.0 1255 


(46a) 


2w„ 1 +// 0.01255 

u h c ~ e 


(46 b) 


Here equation (46a) is an ordinary two-dimensional boundary- 
layer momentum-integral equation for 9 and equation (46b) 
is a relation that evidently must exist among U, bU/bz, 
and 9, when c=0 and w — 0 . 

When the equations are elliptic, no real characteristic 
direction can be found. When only one characteristic direc- 
tion exists, the equations are parabolic and in the hyperbolic 
case two characteristic directions through each point of the 
X 2 -plane are obtained. For the parabolic case then, 


and for the hyperbolic case, 

dz L-\-y'L~ — J(K — J)L 

= U = * W=J) 

dz „ L- y r L^(K=J)L 

(. K-J) 


and the characteristic lines are asymmetric with respect to 
the r-axis. In order to determine whether elliptic, parabolic, 
or hyperbolic equations apply, the magnitude of MN is com- 
puted. Substituting from expressions (30), MN is obtained 
in terms of n : 

6(3n+l)(3n+2) , 1 10 8n + 96 

(5n + 2)(l ln + 7) 55 55(55n 2 +57n+ 14) 

This equation shows MN to be a monotonically increasing 
function of n. For n — 0, 


MN=~ 

and 

Urn MN=^ 

71 — >oo 0 O 

These results indicate that a so-called power-law ] rofile when 
used with £=(l- 2 //$) 2 always results in equations that 
although elliptic are very near to being parabolic. Using 
values (30b), 

MN= 0.936 

which again indicates an elliptic character of the equations. 
It should be remembered, however, that the assumptions for 
G and g were made on the basis of only one set of data; con- 
sequently there is no assurance that the velocity distributions 
existing in the boundary layer will always give the same 
values of MN. In fact, it is generally more likely that they 
will not give the same values of MN, Some indication of the 
variation of MN may already be obtained from figure 8, 



Fioure 8.— Values of parameter MN for (lata of reference 9. 
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/ (re its value was plotted for each Gruschwitz data point. 

points 5, 6, and 7, however, small values of w may have 
re vented an accurate determination of L and as a con- 
equence -MN=0 there. The value of MN in figure 8 
varies within the limits 0.65<MiV<^ 1.2, with the bulk of the 
points indicating that M.V~0.9’5. 

On the basis of the preceding discussion, there is some 
evidence of the equations being parabolic, elliptic, and hyper- 
bolic in the turbulent boundary layer. It is interesting to 
note that generally (as in supersonic and subsonic flow, for 
instance) these hyperbolic and elliptic regions have their 
counterpart in physical phenomena. Thus some essential 
differences might exist in the process of momentum transfer 
between the hyperbolic and elliptic regions. These differ- 
ences cannot now be ascertained because first equations sim- 
ilar to (34) with z ?± 0 would have to be obtained, and there 
is no mention of any irregularities in the behavior of the flow 
in reference 9. When additional experiments are made, 
however, it would seem advisable to study closely these two 
mathematical regions in order to obtain some indication of 
the physical make-up of their differences. 
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